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Using  the  multiple  scattering  Xa  method,  electronic  energy  levels 
have  been  found  self-consis tently  for  the  intermediates  Ni(CO)^,  n=l,2,3 
in  the  formation  of  nickel  t  itracarbonyl  via  the  Ni-CCT  surface  reaction. 
Linear  geometries  have  been  assumed  for  NiCO  and  Ni(C0)?,  and  a  trigonal 
planar  conformation  for  Ni(CO)^,  in  accordance  with  previously  published 
IR  spectra.  The  intermediates  were  assumed  to  be  in  the  gas  phase,  free 
of  surface  interaction.  From  the  energy  level  structure  found,  all  three 
of  the  intermediates  were  determined  to  be  diamagnetic,  as  Ni(CO)^  is 
known  to  be. 
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Introduction 

Self-consistent  field  calculations  using  the  multiple  scattering  Xa 
(MSXa)  technique  are  reported  on  the  electronic  structure  of  the  mono-, 
di-,  tri-,  and  tetracarbonyls  of  nickel.  The  results  allow  the  evolution 
of  the  tetracarbonyl  level  structure  to  be  traced  from  the  point  of  surface 
detachment  of  the  Ni+CO  complex. 

Chemisorption  processes  in  general,  and  catalysis  in  particular,  have 
been  the  subjects  of  intensive  study  in  recent  years.  Of  special  interest 
is  the  role  played  by  d-electrons  in  transition  metal  substrates;  catalytic 
ability  seems,  with  few  exceptions,  to  be  directly  correlated  with  d-band 
vacancies,  filled  by  covalent  bonding  with  adsorbate  molecules.  The  exact 
extent  and  nature  of  d-electron  involvement  in  this  bonding  is  one  of  the 
largely  unresolved  points  of  interest  in  the  study  of  chemisorption.  As  a 
prototype  for  the  understanding  of  nondissociative  chemisorption,  the 
nickel-carbon  monoxide  reaction  has  received  considerable  attention. 

The  bonding  of  the  CO  molecule  to  a  transition  metal  substrate  such  as 
nickel  is  considered  to  be  relatively  simple,  and  therefore  amenable  to 
investigation;  the  CO  on  Ni  system  has  been  extensively  studied  both  exper¬ 
imentally  I-®  and  theoretically Recently,  photoemission  spectroscopy 
of  CO  adsorbed  on  a  nickel  surface  has  determined  the  presence  of  large 
relative  shifts  in  the  CO  levels  upon  adsorption;  the  use  of  angle-resolved 
techniques'^  has  affirmed  that  a  reversal  of  the  two  topmost  occupied  CO 
levels,  the  lx  and  5 a  levels,  takes  place  on  adsorption,  the  5o  level 
becoming  strongly  involved  in  bonding.  As  the  CO  5o  orbital  is  a  carbon 
"lone  pair"  orbital  extending  into  the  metal  surface,  we  would  indeed  expect 
stronger  interaction  than  with  the  CO  lx  orbital  or  the  deeper  oxygen 
"lone  pair"  orbital,  4o,  which  points  away  from  the  surface.  Theoretical 
studies,  concurrently,  have  investigated  the  role  of  nickel  s,  p,  and  d 
electrons  in  bonding  CO,  and  the  extent  to  which  back  donation  of  nickel 
electrons:  into  the  (normally  unoccupied)  CO  2x*  level  takes  place. 

There  have  been  numerous  studies  of  the  Ni(CO)^  molecule2~16,  and,  more 
recently,  calculations  of  the  level  structure  of  atomic  clusters  taken 


to  simulate  CO  bonding  to  a  nickel  surface:  NiCOl?-!^,  Ni^CO^, 

Ni^CO^-^^  an(j  bridge-bonded  N^CO^®-2^ 

Theoretical  approaches  have  included  ab  initio  calculations'^ ^ 
Hartree-Fock-Slater  methods^*22,  and  MSXa  calculations®»^^_^^'^"*-»^^. 

Among  the  calculations  for  Ni(C0)4,  results  as  to  orbital  energies  and 
orbital  character  are  in  agreement  on  general  features;  namely,  the 
highest  occupied  (9 1 2  and  2e)  levels  are  predominantly  of  nickel  3d 
character,  and  the  30  5 o-  and  lr-derived  levels  become  nearly  coincident 
or  reversed.  Calculations  for  chemisorbed  CO  using  cluster  models  also 
predict  a  near-coincidence22  or  reversal-^ 7-21  0f  these  levels. 

There  have  been  several  studies  of  the  linear  NiCO  cluster,  one  of 
the  intermediates  considered  here,  as  a  model  for  chemisorption. 

Bagus  and  Hermann-*-2 » did  Hartree-Fock  calculations  on  NiCO  and  concluded 

•I  T 

that  the  ground  state  had  an  open-shell  (12o  16  )  structure.  Rosen, 

Baerends  and  Ellis,  in  a  subsequent  Hartree-Fock-Slater  calculation^, 
showed  a  closed-shell  configuration  (llo2  ld^) ,  and  recently  Fenske^-*, 
from  a  Hartree-Fock  configuration- interaction  study,  has  also  found  NiCO 
to  be  closed-shell. 

The  primary  motivation  for  the  present  set  of  calculations  was  the 
series  of  experiments  by  R.S.  Mehta2^  investigating  the  effect  of  the 
magnetic  state  of  a  nickel  substrate  on  the  gaseous  Ni(CO)^  formation 
rate.  It  was  found  that  the  activation  energy  /IE  for  the  breaking  away 
of  the  Ni+CO  complex  from  the  surface  is  dependent  on  the  para—  or  ferro¬ 
magnetic  state  of  the  nickel;  for  example,  for  a  Nig  72Cuo  28  a^°y 
/E  *  0.15eV  in  the  paramagnetic  state,  while  £  ~  0.35eV  in  the  ferro¬ 
magnetic  state.  Based  on  these  findings  a  model  for  Ni(CO)^  formation 
was  proposed  in  which  the  Ni+CO  complex  was  inferred  to  be  diamagnetic 
upon  departure  from  the  surface.  Our  calculations  were  undertaken  to  find 
the  magnetic  state  of  the  intermediates  Ni(C0)n,  n=  1,2,3  and  the  evolution 
of  their  electronic  structure  as  additional  carbonyl  groups  are  added. 
Computations  for  CO  and  Ni(f’c/.J  were  included  for  comparison. 


Our  results,  showed  that  all  of  the  series  Ni(CO)n,  n  =  1, 2,3,4  are 
closed-shell  molecules,  anc,  that  most  of  the  principal  level  shifts 
due  to  bonding  occur  by  the  time  the  second  CO  group  has  been  added 
to  the  nickel  atom.  This  :.s  in  accordance  with  the  findings  of  photon 
emission  studies,  which  also  show  large  shifts  in  the  CO  orbitals  for 
a  single  isolated  molecule  bonded  to  Ni. 

Method 

The  multiple  scattering  Xo  method  has  been  extensively  applied  to 
atomic  clusters^,  26  ^  ^  muffin— tin  potential  is  assumed,  constant  be¬ 

tween  atomic  spheres  surrounding  each  nucleus  and  spherically  symmetric 
inside  the  atomic  spheres  andoutside  a  bounding  sphere  enclosing  the 
entire  molecular  cluster.  The  orbital  wave  function  (WF)  is  expanded 
inside  the  atomic  spheres  in  the  partial  wave  representation,  using 
numerically  generated  solutions  to  the  atomic  Schroedinger  equation. 

In  the  intersphere  region  the  expansion  is  in  terms  of  modified  Bessel 
functions.  The  Slater  local  approximation  to  the  exchange  potential, 
Vxch^r)  “  -SapptrJ/Sir]1/3  is  used,  where  a  is  a  numerical  parameter 
of  the  order  of  0.7,  Using  boundary  conditions  of  continuity  of  the 
orbital  WF  and  its  first  derivative  across  sphere  boundaries  the  scattered 
wave  secular  equations  can  be  set  up.  The  generated  expansion  coeffi¬ 
cients  are  then  employed  in  setting  up  the  new  potential,  and  the  process 
is  repeated  to  self-consistency. 

The  method  has  been  widely  used  in  finding  the  electronic  structure 
of  molecular  species  (e.g.,  refs.  8,26)  and  of  clusters  taken  to  simulate 
a  solid^,  in  calculations  of  the  levels  of  NiCCO)^,  those  obtained 
by  MSXa  methods®* are  agreeinent  with  those  found  from  other,  ab 
initio,  methods^*^  with  the  possible  exception  of  finding  a  larger  Ni  4s 
involvement  in  bonding  and  a  lesser  effect  of  2ir  backbonding,  Larsson  and 
Braga^-2  and  Braga,  Larsson,  and  Leite^3  have  shown  convincingly,  however, 
that  the  differences  lie  in  the  interpretation  of  the  resulting  orbital 
WF’s  rather  than  in  the  reliability  of  the  methods  employed. 


All  the  molecules  Included  in  our  study  were  taken  to  be  in  the 

gaseous  phase,  completely  desorbed  from  the  surface.  The  clusters 

consisted,  therefore,  of  a  single  Ni  atom  bonded  to  1,2,3  or  4  CO 

groups.  Calculations  were  carried  to  a  self-consistency  of  ±0.002 

Rydberg  for  all  valence  levels  over  an  energy  range  from  -2.5  Ryd 

to  -0.1  Ryd.  Symmetry  assignments  for  the  intermediates  were  made 

in  accordance  with  the  IR  vibrational  mode  spectra  taken  by  DeKock^®. 

This  gives  a  linear  monocarbonyl  (G^)  and  dicarbonyl  (D®  and  a 

trigonal  planar  (Djj, )  tricarbonyl  as  shown  in  Figure  1.  The  interatomic 

distances  used  for  the  Ni(C0)n,  n  -  1-4  were  those  of  Braga,  Larsson, 

and  Leitex  in  their  calculation  of  the  Ni(C0>4  electron  levels. 

These  distances  were  used  in  the  absence  of  data  on  the  intermediate 

carbonyls.  The  CO  molecule  equilibrium  distance  was  used  for  the  isolated 
29 

CO  calculation.  The  radii  of  the  atomic  sphere  regions  were  chosen 
to  have  a  20%  overlap  to  improve  accuracy;  previous  studies-^* 31  ^ave 
concluded  that  such  an  overlapping-spheres  model  provides  an  approximate 
correction  for  non-muffin  tin  effects  by  reducing  the  volume  of  the 
constant-potential  intersphere  region.  The  carbon  and  oxygen  atomic 
sphere  radii  were  determined  from  the  ratio  of  covalent  radii  given 
by  Slater-^.  Atomic  cores  were  not  frozen. 

Results  and  Discussion 

The  calculated  nickel  monocarbonyl  valence  levels  are  shown  in 
Figure  2  in  comparison  with  the  CO  levels.  The  isolated  CO  calculation 
yielded  level  orderings  and  characters  in  agreement  with  the  known  CO 
level  structure.  Dashed  lines  connect  the  monocarbonyl  levels  to  the 
CO  levels  from  which  they  are  principally  derived.  To  differentiate  the 
CO  levels  from  those  of  NiCO,  which  have  the  same  symmetry  labelling, 

NiCO  levels  will  be  indicated  in  the  following  discussion  by  a  tilde, 
e.g.  99.  The  lowest  NiCO  valence  levels,  89  and  99,  are  the  shifted  CO 
3a  and  4a  Levels,  respectively.  The  CO  In  level  shows  a  shift  of  only 
.03  Ryd,  or  0.4eV,  upon  becoming  3ft,  while  the  5o  level  shifts  downward 
by  3.4eV  to  become  109,  the  i  incipal  NiCO  bonding  level.  This  consti¬ 
tutes  a  reversal  of  0.9eV  between  the  109  a^d  3fl  levels;  angle-resolved 
photoemission  measurements  give  a  0.5eV  reversal  for  CO  on  Ni(100)^  and 


a  1.6eV  reversal  on  Ni(lll)^,  by  way  of  comparison.  As  we  will  see, 
the  addition  of  successive  CO  groups  to  Ni  widens  this  reversal  separa¬ 
tion.  The  105  level  is  an  admixture  of  the  5o  level,  the  carbon  "lone 
pair"  orbital,  with  Ni  s,  pz,  and  dz2  orbitals  to  a  virtually  equal 
extent  (31%  s,  31%  pz,  and  38%  dz2). 

The  highest  occupied  NiCO  orbitals  are  the  115  and  16  .respectively, 
both  nickel  d-nonbonding,  a  feature  that  persists  throughout  the  series 
of  intermediates.  The  4tt  level  is  nickel  d  (90%)  with  a  CO  u  (10%) 
admixture.  There  are  no  partially  filled  orbitals  in  NiCO;  thus  it  has 
no  static  magnetic  moment  according  to  our  calculation.  Had  the 
(empty)  level  been  more  nearly  degenerate  with  the  116  and  l5  levels, 
one  might  have  expected  exchange  splittings  to  give  partially  filled 
levels  and  thus  a  paramagnetic  molecule.  The  gap  between  115  and  5  if, 
however,  is  large,  about  2eV,  so  the  molecule  is  diamagnetic. 

Dicarbonyl  valence  levels  are  compared  to  those  of  CO  in  Figure  3. 

The  CO  3o  and  4o  levels  are  split  into  the  degenerate  (to  the  accuracy 
of  our  calculation)  pair  6o+,  5cr±  and  the  close- lying  pair  7o+,  6o£, 
respectively.  The  CO  lif  level  participates  in  the  formation  of  three 
occupied  valence  pi- levels,  the  3ttu,  litg,  and  2irg.  The  3itu  orbital 
is  virtually  completely  CO-like,  showing  only  about  1%  Ni  character. 

The  lffg  and  2wg  orbitals  represent  combinations  of  the  CO  lir  and  2tt* 
vavef unctions.  liTg  is  a  pure  oxygen  p-type  orbital,  with  no  nickel 
or  carbon  character,  while  2Tig  has  no  admixture  of  oxygen  and  combines 
carbon  p-character  with  Ni  dKZ)yZ.  Thus  the  liTg  and  2itg  orbitals  provide 
the  first  direct  evidence  of  occupation  of  the  CO  2tt*  orbital  in  the 
Ni(CO)n  series. 

The  CO  5o  level  participates  in  forming  the  80g+,  7ou+  orbitals, 
which  again  lie  below  the  lir-derived  levels.  The  reversal  splitting 
between  1ttu,  the  most  "CO  lTr-like"  orbital  and  the  averaged  8og+,  7ou+ 
pair  is  now  2.9eV,  a  2eV  increase  over  the  NiCO  analogous  separation. 

The  topmost  occupied  orbitals  in  Ni(C0)2  are  the  9ag+  and  16^  nonbond¬ 
ing  levels,  as  expected.  As  these  are  filled,  the  dicarbonyl,  like  the 
monocarbonyl,  is  diamagnetic. 


The  final  molecule  studied  was  the  tricarbonyl t  Figure  4t  The  same 
general  features  are  evident  here  as.  in  Ni(CQ)2  and  NiCO,  The  3o  and  4o 
levels  of  CO  split  into  e*,  a^’  pairs;  the  set  of  levels  7a2%3a2", 

8e’  and  7e" , clustered  around  -0.81  Ryd,  comprise  the  CO  ln-evolved  levels 
(while  2n*  involvement  can  be  assumed  here, given  the  backbonding  present 
in  Ni(C0)2  and  Ni(CO)^  it  cannot  be  unambiguously  demonstrated  within 
the  limits  of  our  calculation).  The  5o  level,  split  into  8aj’  and  7e’, 
is  an  average  of  2.6eV  below  the  cluster  of  ln-related  levels.  Nickel 
d-character  is  distributed  principally  on  the  set  of  levels  9e'',  2e", 
and  Sa^',  the  highest  occupied  orbitals.  Again  in  NiCCO)^  all  occupied 
orbitals  are  completely  filled. 

It  is  of  interest  that  our  results  confirm  to  some  extent  DeKock's 
.i^signment  of  molecular  geometries.  Given  that  NiCO  and  Ni(CO)^  are 
closed-shell,  the  fact  that  and  03^  symmetries  lead  to  closed-shell 
conf igu'ations  for  Ni(C0>2  and  Ni(C0)3  tends  to  substantiate  these  symmetry 
assignmei.t? . 

Figure  5  depicts  the  evolution  of  the  level  structure  from  CO  to 
Ni(CO)^,  showing  only  the  principal  CO-derived  valence  levels.  The 
calculated  NiCCO^  level  structure  is  in  close  agreement  with  that  found, 
for  instance,  by  Braga,  Larsson  and  Leite^  in  their  MSXa  calculation;  we 
will  not  elaborate  here  except  for  purposes  of  comparison  with  an  inter¬ 
mediate  species.  The  energy  difference  between  the  CO  lu-derived  levels 
8t2>  and  le)  and  the  principal  5o-derived  levels  (7t2»8a^)  is  2.8eV 
on  the  averagein  Ni(CO)^.  (One  should  not  put  undue  emphasis  on  these 
estimates  of  the  reversal  shifts,  as  they  are  simple  averages  over  the 
groups  of  levels  containing  significant  In-  and  5o-character .  They 
serve,  however,  to  illustrate  the  trend  as  CO  groups  are  added.)  The 
following  general  conclusions  can  be  drawn  from  comparison  of  the  Ni(C0)n, 
n  =  1-4  series:  1)  the  low-lying  CO  3o-  and  4c-derived  valence  levels 
shift  very  little,  as  expected,  during  the  formation  of  Ni(CO)^  (the  3o 
level,  for  instance,  going  f ro  1  -2,16  Ryd  in  CO  to  -2.21  Ryd  in  Ni(C0)2» 
and  -2.14  Ryd  in  Ni(C0)4);  2)  the  nickel  nonbonding  levels  shift  slightly 


toward  lower  energies  as.  the  number  of  CO  groups  attached  increases 
(see  Table  II);  3)  most  of  the  shifts  associated  with  the  CO  lrc  and 
5o  levels  have  occurred  by  the  time  the  second  carbonyl  group  has  been 
added,  the  reversal  shift  growing  from  0,9eV  in  NiCO  to  nearly  3eV  in 
NKCO)^. 


Summary 

It  has  been  cc.ifirmed  that  the  isolated  Ni(CO)n  molecule  is  diamag¬ 
netic  for  any  value  of  ri  from  1  to  4,  Thus  the  Ni+CO  complex  is  expected 
to  be  nonmagnetic  upon  detachment  from  the  Ni  substrate  during  the  Ni¬ 
CO  surface  reaction.  Reversal  of  the  5  and  1  CO-derived  levels  takes 
place  even  with  only  one  attached  CO  group  for  the  interatomic  distances 
used  in  the  present  calculations. 
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Table  I 


C  a 1 culational  Parameters 

o  •  factors 

carbon:,  0.75928 
oxygen:  0.74447 
nickel:  0.70896 
intersphere  (for  CO) :  0 

intersphere  (for  Ni(C0)n 

Sphere  Radii  (in  a.u.) 

CO:  carbon  1.378 
oxygen  1.181 
or ter  sphere  2.346 

Ni (C0)n: 

carbon  1.380 
oxygen  1.182 
nickel  2.749 

outer  sphere  for  NiCO  4.753 

outer  sphere  for  Ni(CO)2_4  6.757 


Interatomic  Distances  (in  a.u.) 

R  for  CO:  2.132 
co 

Rnic  for  Ni(C0)n:  3.440 
Rco  for  Ni(C0)n:  2.135 


.75039 
):  0.72376 


Table  II 


Final 

Energy  Levels  for 

CO  and  for  Ni(CO)^,  n=l-4 

(Rydbergs) 

Levels  above  dashed  lines  (_ 

_)  are  unoccupied. 

CO 

NiCO 

Ni(C0) 0 

+  L  . 

— 

_  _-0.199 

12o  -0.223 

5  a 

-0,686 

5tt  -0.262 

4ttu  -0.290 

llT 

-0,899 

llo  -0,412 

3h£  _  _-0.380 

4o 

-1,088 

16  -0.415 

9ag  -0.543 

3o 

-2.088 

4ir  -0.434 

16g  -0.584 

2a(Cls2) 

-20.161 

3  it  -0.871 

2rg  -0.637 

lo(01s2) 

-37.798 

10 a  -0.938 

lirg  -0.659 

total 

9o  -1,194 

3nu  -0.811 

energy 

—224.751 

8a  -2.227 

7ou+  -0.972 

7a,2ir(Ni3p*’) 

8ag+  -1.089 

-4.882 

o 

6ou+  -1,162 

6a(Ni3s  )  -7.593 

o 

7ag+  -1.203 

5a(Clsz)  -20.383 

9 

5agV  -2,208 

4a(01sz)  -37.742 

& 

6ag+  -2.208 

3a, 1  (Ni2p  ) 

-61.592 

4au+,2au  (Ni3p^) 
-5.060 

2a(Ni2s  )  -70.495 

2 

5ag+(Ni3s2)-7.772 

la (Nils  ) -594. 322 

3au+,4ag+(Cls2) 

total 

-20.476 

energy  -3236.255 

2au+,3og+(01s2) 

-37.613 

iau+,luu(Ni2p6) 

-61.784 

2og+(Ni2s2) 

-70.652 

log+(Nils2) 

-596.500 


total 

energy  -3457.366 


Tabic  II.  continued 


Ni (CO) ^ 


3e" 

-0,047 

lie' 

-0,058 

3a2 

-0,085 

4a'2 

▼0.137 

lOaj 

-0.155 

10e' 

-0.176 

2a2 

-0,314 

9e' 

-0.464 

9a| 

-0,496 

2e" 

-0.508 

3a2 

-0.800 

lai 

-0.811 

le" 

-0.811 

8c' 

-0.816 

7e' 

-0.972 

8a' 

-1.023 

6e' 

-1,55  — 

7ai 

-1,164 

5e' 

-2.175 

6a' 

-2.175 

1 

2a’^4e' 

(Ni3p6) 

t.4.961 

5a|(N13s2)  -7.673 

4al,3e' (Cls2) 

-20,379 

3a’,2e'(01s2) 

▼37.716 

la",le'(N12p6) 

-61.690 

2a^  (Ni2s2)-70. 559 

la' (Nils2) 

-596,399 

total 

energy  -3680.925 


Ni(C0)4 


luz 

-0.065 

-0.102 

9al 

-0.159 

3c 

-0.168 

iot2 

/-0.196 

9t2 

-0.490 

2e 

-0.558 

lfc1 

-0.781 

<N 

U 

00 

-0,801 

lc 

-0.807 

7t2 

-0.966 

8al 

-1.021 

-1.129 

-1.151 

5t2 

-2.137 

631  6 
4t2iNi3p  ) 

5u^(Ni3s2) 

4a,3t2(Cls 

-2.137 

-5.013 

-7,727 

2) 

-20.325 

3a,2t,(01s2) 

z  -37,650 


lt2(Ni2p6)-61.756 

2a],(Ni2s2)-70.626 

la  (Nils2) 

-596.448 

enerfy  -3904.945 


a a  (co) 


3  He! 

3c-  "  //  ^  2 

—  10  C. 

IOC 

- Za'z 


Hc' 


- v«; 

Ze'1 


cccn  piedi 


r 7  t 


<o<z»,  5e' 


Figure  4 


Rydberg 


